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Abstract 

We discuss the role of infrared renormalons in power corrections to hadronic cross 
sections, including their universality. We show how perturbative renormalon structure 
arises near kinematic boundaries in the thrust distribution, the Drell-Yan cross section 
and radiative B decays. The leading infrared renormalon in each case is associated with 
jet evolution. Demanding that the combination of perturbative and nonperturbative 
contributions be well-defined, we infer the form of the leading power corrections to each 
cross section. This is at the level 1/Q for the thrust distribution and Drell-Yan process, 
and 1/m? for B decay. We discuss the universality of 1/Q corrections between Drell- 
Yan and thrust, and conclude that it is approximate, due to contributions from multijet 
configurations that may differ in the two cases. 
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1. Introduction. Our current understanding of hard processes in QCD is based on factor- 
ization (the "improved the parton model"), which allows us to calculate cross sections in the 
high energy limit Q 2 — > oo (with Q 2 a large kinematic scale) as a convolution of nonpertur- 
bative parton distribution functions and short distance partonic subprocesses. The latter are 
calculable in perturbative QCD as a series in a s (Q 2 ). Although this procedure is adequate 
for a variety of processes, there is a wealth of data for which low-order perturbation theory 
alone is not enough, most notably the transverse momentum distributions of W ± and Z° in 
the Drell-Yan process 0, 0, [|, || and direct photon and diphoton production || || , as well as 
event shapes in e + e~-annihilation J7|, ^J. As we approach the edge of the phase space in these 
processes, we observe a disagreement with next-to-leading perturbative QCD calculations that 
can be attributed to nonperturbative (or hadronization) effects, our understanding of which is 
still incomplete. Remarkably enough, in some cases one may achieve a better agreement with 
experimental data by including additional factors into parton model predictions "by hand". 
These factors take into account nonperturbative effects in the form of power corrections in 
l/Q. In this paper we would like to describe a simple method |§ which allows us to understand 
the appearence of additional nonperturbative factors in resummed cross sections. It is based 
on the analysis of ambiguities of the perturbative series in QCD associated with the so-called 



infrared renormalons HO, IIT 



In the special case of processes (like the total cross section for e + e~ — > hadrons, deep- 
inelastic scattering, etc.), which admit the operator product expansion (OPE), nonperturba- 
tive power ("high twist") corrections to a physical distribution a(Q) can be parameterized by 
dimensional nonperturbative scales a n as 

oo 

°(Q) = v P cAQ) + E %■ (!) 

n=no ^ 

Here the first term corresponds to the perturbative contribution cr per t = J2k c k(o; s (Q 2 )) k and the 
leading power corrections appear as ~ Q~ n ° . While the complete sum in (|I|) defines a physical 
quantity a(Q), individual terms are not well defined, either due to ambiguities in the definition 
of nonperturbative quantities a n ~ Aq CD or due to divergences of <7 per t related to the factorial 
growth of the coefficients Ck ~ P^kl (with f3 Q — 11 — 2/3n/) at higher orders of perturbative 
QCD. As a result, to give a meaning to the perturbative (as well as nonperturbative) contribu- 
tion to a(Q) one has to specify the way in which we "regularize" divergences of cr pert . Different 
prescriptions lead to different results for cr pcrt which eliminate ambiguities at the level of power 
corrections AQ CD /<5 n . It is important for us that ambiguities of the perturbative contribution 
crpert are compensated in ([!]) by ambiguities in the definition of nonperturbative scales a n . Using 
this property we can explore the ambiguities of <7 pert to identify the structure of nonperturba- 
tive corrections, and in particular the level n at which the leading power corrections appear in 
([I]). A justification for this scheme is partly in our expectation that there is a unique theory of 
QCD involving perturbative and nonperturbative effects, and also that for processes in which 
the OPE is applicable and no is known explicitly, the analysis of perturbative ambiguities leads 
to the same result |I3J . 

Let us apply the same approach to hadronic processes to which the OPE is not applicable. 
To be concrete, we will consider the following three different processes: the thrust distribution 
da/dT in e + e~ annihilation, the energy distribution of lepton pairs in the Drell-Yan process 
and the differential rate of the radiative decay of the B meson. One of the reasons why the 
OPE does not apply directly to these processes is that they involve two scales Q 2 and q 2 with 
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Q 2 ^> q 2 ^> Aqcd- (We shall see below, however, an important role for the OPE in heavy quark 
effective theory |12|.) As a result, the short-distance partonic subprocess includes large double- 
logarithmic (Sudakov) corrections (a s (Q 2 ) In 2 (Q 2 / q 2 )) n , due to the enhancement of soft gluons 
with momenta k ~ q. Our strategy will be the following |J. We first apply factorization, and 
perform resummation of large perturbative logarithmic corrections. The resulting resummed 
expressions involve integrals of the running constant at very low scales, given by soft gluon 
transverse momenta. This makes the perturbative series sensitive to the way in which we 
treat singularities of the coupling constant. Once we identify the ambiguities of the Sudakov 
resummed perturbative expansions, we add nonperturbative corrections in order to compensate 
them and make the final expressions well-defined. Although the analysis of perturbation theory 
cannot give us the absolute normalization of nonperturbative corrections, we may use it to 
parameterize nonperturbative effects by introducing new scales with dimensions fixed by the 
positions of infrared renormalons. 

2. Sudakov resummation. For e + e~ annihilation, the thrust is defined in terms of the 
momenta kj of final state particles as 

T = max _ ?lh^3 . (2) 
n £;l%l 

It takes a maximum value T = 1 for two infinitely narrow jets. For 1 — T small the final state 
consists of two jets, with invariant masses k 2 and k 2 , respectively, and with energy Q in the 
center-of-mass frame, Q 2 ^> k 2 , k 2 . Then, in the T — > 1 limit the thrust is given by 0] 

P P 

IX fx, , „ 

and for the thrust distribution one can derive the factorized form [[?] 

T - l? k i: dp MQ^)MQVk>) s(i -T- 1 - 1) , (4) 

where <7 to t is the total cross section of e + e~ — > hadrons. Here, the jet distribution J q (Q 2 ,k 2 ) 
describes the probability for a quark with energy Q and invariant mass k 2 -C Q 2 to create a jet 
of collinear particles accompanied by soft gluon radiation ]13|, In the Born approximation, 
= S(k 2 ), while the emission of soft and collinear particles by the "parent" quark leads to 
large logarithmic corrections a™/k 2 ln 2n_1 Q 2 /k 2 to J q . They can be resummed using evolu- 
tion equation techniques, and the final expression for J q can be represented to next-to-leading 
logarithmic accuracy in the form of a Laplace transform 0, |14], |^], 



rdk 2 e-^^ 2 J q (Q 2 ,k 2 )=ex P \-f -(l-e— ) f cusp {a s {k 2 t )) + iMuQ 2 )) 

Jo [ Jo U v ' [Ju 2 Q 2 kt 

(5) 

An identical expression applies to the antiquark jet distribution. Here, v is the transform 
parameter, and r cusp (a s ) and 7(a s ) are anomalous dimensions, known to the lowest orders of 
perturbative QCD. We notice that the jet distribution @ satisfies the normalization condition 

/■oo 

/ dk 2 J q (Q 2 ,k 2 ) = 1, (6) 
Jo 
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which guarantees that / dT^ = a tot . This means, that in agreement with the OPE for 
Ctot(Q); large perturbative (as well as nonperturbative) corrections to the differential cross- 
section da/dT described by (f|) and (|5|) cancel in the total cross-section. Using (£|) we find the 
following expression for the thrust distribution 



-iz(l-T) 



1 rrr, 1 da 

dT — — e 
a dT 



-u(l-T) 



poo 

J dk 2 e- vk *l Q2 UQ\k 7 



(7) 



Together with (^j) this relation resums all large perturbative corrections which describe the 
effect of Sudakov suppression of the end-point region T ~ 1 in the thrust distribution da/dT. 
We notice, however, that in (|5|) the integration over the transverse momenta of soft gluons 
is potentially divergent due to singularities of the coupling constant a s (fef), and it is this 
property of the Sudakov resummed perturbative expression which we are going to use to identify 
nonperturbative effects. 

As a second example, we consider the distribution of the lepton pairs in the Drell-Yan 
process h\ + h 2 — * + X with respect to their invariant mass Q 2 at the edge of phase space 
as Q 2 approaches the invariant energy of the incoming hadrons ^/s. As r = Q 2 /s — > 1, the 
final state in the partonic subprocess is dominated by soft gluons with total energy ~ (1 —t)Q, 
which lead to the appearence of large Sudakov perturbative corrections to the short- distance 
partonic cross-section |TB|, [T4], [IJ|. The proper quantities to study in this case are the moments 
of the differential cross-section da DY /dr, normalized to the moments of a structure function of 
deep- inelastic scattering F(x), 

, da^ r „ 



2 

n i 



a,, 



dTT T - 



dr 



dxx n F(x) 



Each of these quantities may be factorized into the product of moments of a parton distribution 
function and moments of a short- distance partonic cross-section. The parton distributions 
cancel in the ratio, A n , which makes it possible to evaluate A n perturbatively. For large n, 
the ratio A n has large Sudakov perturbative corrections (a s ln n) k , whose origin can be traced 
back to the asymptotics of da DY /dr and F(x) as r — > 1 and x — > 1, respectively. In the Drell- 
Yan process, large corrections to da DY /dr are caused by soft gluons emitted from the initial 
state by the quark and antiquark before they annihilate. At the same time, in deep-inelastic 
scattering large corrections to F(x) as x — > 1 arise from soft gluon emission from the incoming 
quark, as well as from the decay of the scattered quark with energy ~ Q and small invariant 
mass k 2 = Q 2 (l — x)/x into a jet of collinear and soft particles in the final state. We recognize 
that the latter subprocess is described by the same jet distribution J q (Q 2 ,k 2 ) which entered 
into the analysis of the thrust distribution. Moreover, in the ratio of the moments A n the 
contributions of soft gluons emitted by quarks in the initial state in the Drell-Yan process and 
deep-inelastic scattering cancel, and the only contribution which survives is that of the jet 
distribution associated with the outgoing quark in deep- inelastic scattering. As a result, for 
large n, A n (Q) is given by 



A n (Q) 



H(a s (Q 2 )) 
H(a s (Q 2 )) 



Q 2 [ dxx n JjQ 2 ,Q 2 (l-x)/x) 
Jo 

dk 2 e~ nk2 / Q2 UQ 2 ,k 2 



-2 



1 -2 



0{l/n) 



(9) 



where H = 1 + J2k hh a s(Q 2 ) comes from the contributions of hard virtual gluons with momenta 
k ~ Q, and where the jet distribution has been defined in (^|). Comparing this expression 
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with fl7|) we notice that the resummed expressions for different physical quantities, A n (Q) and 
e -n(i-T)\ fa e same form. 

As a last example we consider the differential rate for the inclusive radiative decay B — > 7X3 
in the end-point region of the photon spectrum, in the limit when the mass of the s quark is 
neglected. In the rest frame of the B meson, we define a scaling variable x as the ratio of 
the photon energy to the mass of the b quark, x = 2E 1 /m. That is, we are interested in the 
inclusive distribution dT/dx with x ~ 1. In the heavy quark limit, m — * 00, we may factorize 
the decay distribution into a convolution of the b quark distribution function in the B meson 
and a short distance partonic subprocess b — > sj. The latter gets large perturbative corrections 
~ [ct"/(l — x)] ln 2n_1 (l — x) in the end-point region, which originate from the propagation 
of the s quark into the final state with a large energy ~ m/2 and a small invariant mass, 
k 2 = m 2 (l — x). The s quark creates a jet of collinear particles accompanied by soft radiation, 
which is described by the same jet distribution J q (m 2 ,m 2 (l — x)) that we encountered above. 
Introducing the moments of the differential rate dT/dx as 

/■Zmax 1 dF 

M n (B^ 1 X s ) = j Q dxx n f —, (10) 
and performing the Sudakov resummation of large perturbative corrections to A4 n , for large 



moments n we may represent the result in the following factorized form [16] 



M n = H(a s (m 2 )) J n S n + 0(l/n) , (11) 

where fff> is a nonperturbative parameter describing the distribution of the b quark in the B 
meson, and H comes from hard gluon corrections to the effective vertex b — ► S7. It is similar 
to the analogous factor in @. J n takes into account all effects of the s quark fragmentation 
into a jet and is equal to 
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./„ nr l I dxx n J q (m 2 , m 2 (l - x)) = f°° dk 2 e~ nk2 /m " JJm 2 , k 2 ) + 0(l/n) , (12) 
Jo 



with the jet distribution given by fl5|). There is an additional factor S n in (|1T[) which takes into 
account the nonleading contribution of the soft gluons emitted by the incoming b quark. It is 
defined as 

^ = exp|-| o 1 ^(l-e-™)r(« s ( M 2 m 2 ))| , (13) 

with r(a s ) an anomalous dimension flIE| . Comparing ( p!3|) and (|5|), we notice that the contri- 
bution of S n to the moments J\4 n is subleading with respect to that of jet subprocess, J n , so 
that it can be reabsorbed into a redefinition of the anomalous dimensions r cusp and 7 to higher 
orders of perturbation theory. 

Thus, for three processes considered above the Sudakov resummation of large perturbative 
corrections gives similar expressions, (0), (|9|) and (|TT|). The reason for this universality is that 
in all cases large perturbative corrections came from the same jet subprocess which can be 
described by the distribution function J q defined in (|). 

3. Power corrections from infrared renormalons. Expressions (0), (0) and ( |TT| ) were 
found after resummation of large Sudakov logarithmic corrections to all orders of perturbation 
theory. This does not guarantee, however, that the resummed perturbative series are conver- 
gent. Indeed, it is clear from (H) that the integral over the transverse momenta of soft gluons 
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is divergent at small values of fcf. In the perturbative expansion of the exponent in the r.h.s. 
of dD in powers of a s (Q 2 ) these divergences manifest themselves in factorial growth of the 
coefficients ~ fi^kl. This tells us that, taken alone, resummed perturbative predictions (0), (|9|) 
and (|TT|) are not well defined. Since this problem arises from small momenta region of soft 
gluons, we do not expect perturbative QCD alone to suffice. Rather, perturbative expressions 
should be supplemented by nonperturbative corrections, and it is their sum that should give 
unique predictions for physical quantities. Although we do not know in general how to evaluate 
nonperturbative effects from first principles, we may explore the ambiguities of perturbative 
expressions like ([5]) to indentify potential sources of nonperturbative corrections. 

Let us identify the leading infrared renormalon singularity in the expression for the jet 
distribution (|5J). It is defined by the first term in the exponent, in which we interchange the kt 
and u integrals and expand the result for small transverse momentum to get |9j 



exp 



2v 

~QJo 



dh T cusp (a s (k 2 ))[l + 0(k t /Q)} 



(14) 



We conclude that the resummed perturbative expression for the jet distribution (^) has its 
leading infrared renormalon singularity at the level of a 1/Q power correction, 



r dk 2 e- vk2 l Q2 JJQ 2 ,k 2 
Jo 



nonpert 



exp (-vA/Q + 0(1/Q 2 



(15) 



where the dimensional parameter A includes both perturbative and nonperturbative contribu- 
tions. Although the latter may depend on the process in which the jet distribution enters, the 
former have a universal form (~ 2 dk t T cusp (a s (k 2 )), which follows from flTil) . Applying (15 ) 
to (0) and ([|) we estimate the magnitudes of nonperturbative effects in the thrust distribution 
and in the Drell-Yan process as 



-Ki-T) 



~ e 



-2vA T /Q 



nonpert 



A n (Q) 



IQ 



nonpert 



(16) 



This analysis of perturbative ambiguities, and their relation to nonlocal operators, suggests |Q 
that the contribution of the leading infrared renormalon to the parameters A^ and A DY is the 
sameQ, although it does not allow us to draw definitive conclusions on the magnitudes of these 
parameters. 

One may try to apply the same result to the moments of the B meson differential rate, 
eq.(|TTJ). There is, however, an important difference from the previous cases. For B decay, the 
energy scale Q, which is a kinematic invariant for the cross sections (center of mass energy 
for thrust and lepton mass in the Drell-Yan process), should be replaced by the b quark mass 
m. The 1/m contribution is given by an operator expectation value in heavy quark effective 
theory that vanishes by the heavy quark's equation of motion. This result may also be seen 



directly in perturbation theory. It turns out QT9| that the definition of m is also ambiguous 
due to the presence of an ultraviolet renormalon in the perturbative series for m. Therefore, 
analyzing perturbative ambiguities to the moments J\4 n of the differential rate, we have to take 



into account contributions of infrared renormalons to the jet distribution fll5|) and that of the 
ultraviolet renormalon to the b quark mass, m — > m + Sm. The net effect of the latter ambiguity 



1 This result has been recently confirmed in 0, ^8) 
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can be reduced to a renormalization of the scaling variable x — > x(l — 5m/m) in (|T0|) and the 
renormalon contributions to the moment M. n can be represented in the large m limit as 



M n 



~ exp(— nA/m — n5m/m) (17) 

nonpert 



with parameter A defined in (|15|). The absence of 1/m corrections implies that their sum 
vanishes, A+5m = 0, and perturbative ambiguities will first appear in M. n at the level 0{l/m 2 ) 
power corrections. 

We stress that although the resummed perturbative expressions for the thrust distribution 
and in the B meson decay, eqs. (|^) and (|TTD, look very similar after identification of the 
energy Q with heavy quark mass m, their properties with respect to renormalon ambiguities 
are completely different, due to ambiguities in the definition of heavy quark mass m. However, 
keeping in mind the situation with B decay, we may interpret nonperturbative 1/Q corrections 
to the thrust distribution as a renormalization of the scaling variable 1 — T, or equivalently 
the invariant masses k 2 and k 2 of the jets. This implies, in particular, that the emission of 
nonperturbative soft gluons by the quark and antiquark makes the jets in the final state wider, 
(k 2 + k 2 )/Q 2 > At/Q, and the phase space for the thrust becomes smaller, T < 1 — At/Q- 
Indeed, using the expression (0) for the thrust distribution, and expanding it in powers of z/, one 
can estimate the nonperturbative correction to the average thrust from two-jet configurations, 



(i-r> 



2 rl ct A T /Q. (U 

nonpert 



In contrast to the thrust distribution ([/]) at v ^> 1, the average thrust (1 —T) gets leading 
contributions not only from 2-jet events but from multi-jet events as well. The latter should 
be analyzed separately, and there is no reason to expect that their contribution to the average 
thrust will be subleading by a power with respect to that in (0). Indeed, an analysis of 
renormalon ambiguities to perturbative series corresponding to multi-jet final state shows |^TJ 
that leading nonperturbative corrections to (1 — T) appear at the same level 1/Q as in the 
2-jet events, although suppressed in general by powers of a s {Q). Thus, the total leading 
nonperturbative correction to the average thrust has the same form (|18D but with the scale 
At modified by the contributions of multi-jet events. Contributions of infrared renormalons in 
1/Q— power corrections to (1 — T) associated with multi-jet events prevent an exact equality 
between perturbative ambiguities in the average thrust and in the Drell-Yan process. 

4. Conclusions and perspective. The analysis of ambiguities associated with infrared renor- 
malons allows us to gain insight into the structure of nonperturbative corrections in hadronic 
processes to which the OPE is not applicable. Performing Sudakov resummation of large per- 
turbative corrections, we found that in all three cases, the thrust distribution, the Drell-Yan 
process and the differential rate of the B meson decay, the contribution of the leading infrared 
renormalon has a universal form, which suggests that nonperturbative effects may appear at 
the level of 1/Q (or 1/m for the B meson decay) power corrections. Moreover, the 1/Q— power 
corrections to the thrust distribution and to the Drell-Yan process can be parameterized by 
introducing two new nonperturbative scales At and A DY , which in an analogy with the local 
condensates in the QCD sum rules can be related to vacuum expectation value of nonlocal 
Wilson line operators |9|, |21|. An important difference occurs for B meson decay, where the 
OPE in heavy quark effective theory ensures that infrared and ultraviolet renormalons to the 
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1/m power corrections cancel each other, indicating that leading nonperturbative effects to the 
differential rate occur at the level of 1/m 2 — power corrections. 

The above considerations can be easily generalized to other processes which require resum- 
mation of large perturbative corrections. Returning to the example mentioned at the very 
beginning, we consider the transverse momentum distribution || of lepton pairs in the Drell- 
Yan process d 2 a / dq 2 dQ 2 . In contrast with the previous case, we do not require Q 2 ~ s but 
approach the edge of the phase space g ( 2 C Q 2 . In this kinematic region, the small transverse 
momentum of the lepton pair is compensated by soft gluons emitted by the quark and anti- 
quark in the partonic subprocess. These soft gluons give rise to large perturbative corrections 
«g \n 2k (Q 2 / q 2 ) , which can be resummed to all orders of perturbation theory into a Sudakov 
form factor |IJ. As above, the resulting resummed expressions for the transerse momentum 
distribution suffers from infrared renormalon ambiguities which give rise to nonperturbative 
corrections. The leading nonperturbative corrections have a simple physical meaning |§. They 
introduce an additional broading into the transverse momentum k t — distribution of incoming 
quarks with a gaussian weight. In this fashion, we arrive at the form originally proposed by 
Collins and Soper for nonperturbative contributions to the Drell-Yan (^-distribution |T| , 

2 i Q 



Here, the width a depends on the invariant energy of the quark and antiquark in the par- 
tonic subprocess, s = X1X2S = Q 2 , and the numerical values of new nonperturbative scales 
(gx,g2,Qo) have been estimated M from comparison with experiment as g\ = 0.15 (GeV) 2 , 
Q2 = 0.40 (GeV) 2 and Qo = 2 GeV. As we have seen, the infrared renormalons contributing to 
1 / Q-power corrections have universal structure in different processes. The same is true for the 
infrared renormalon arising in the Sudakov resummation for the q t — distribution. One can also 
show, for example, that resummation of soft gluon effects for direct photon production at small 
xf = 2q t /y/s and for diphoton production at small total q t of the photons, leads to expres- 
sions in which the contribution of the leading infrared renormalon gives rise to nonperturbative 
corrections similar to those in (|I9). 
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